Gap symmetries from the neighbor coupling in square-lattice superconductors 
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The gap symmetries of superconductivity are studied in this work. It is found that the gap 
symmetries are simply determined by the 4-fold rotational symmetries of the coupling potential on 
neighbor sites. A local on-site coupling potential results in the on-site pairing with the conventional 
s-wave symmetry, but a coupling potential between the nearest neighbors or the next-nearest neigh- 
bors results in the pairing on neighbor sites with the s~ , d x 2_ y 2, d xy , or s x 2 y 2 gap symmetries. 
It is proved that both isotropic and anisotropic gap functions are allowed by the 4-fold rotational 
symmetries of the coupling potential. Finally a numerical computation is performed to demonstrate 
the gap symmetries. This neighbor coupling provides a unified picture for the gap functions of the 
conventional and the high Tc superconductivity. 

PACS numbers: 74.20.-z 
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INTRODUCTION 



A few decades of study on high-Tc superconductor con- 
firmed the d-wave symmetry of the gap function at least 
in cuprate superconductors, relative to the conventional 
s-wave BCS gap symmetry P, The observation of the 
half quantum of magnetic flux provides a direct detec- 
tion to the d-wave gap symmetry in cuprates[2]. The 
concept of <i-wave order, however, is quite misleading to 
the physics of superconductivity. In general, it is inferred 
from the antisymmetry requirement of the wave function 
of a Cooper pair. Since electrons are fermions the wave 
function of a pair of electrons must be antisymmetric un- 
der their exchange. When the spin part of the wave func- 
tion is antisymmetric the real space part should be sym- 
metric. Then it was thought that the angular momentum 
of a Cooper pair can only be even, i.e., L — 0,2,4,.... 
These quantum numbers led to the description of s, d, g- 
wave order parameters. This picture is a good approx- 
imation for electrons on a atom, where the rotational 
isotropy of space allows the angular momentum to be a 
good quantum number. When the coupling potential be- 
tween electrons is attractive on one atom it can be indeed 
expected to observe these gap symmetries. Nevertheless, 
for high-Tc superconductors, the coupling is no longer 
on-site attractive, but is attractive only on neighbor sites. 
In this case the space isotropy is broken by the crystal lat- 
tice and only the point group symmetries remain. There- 
fore, the concept of d-wave symmetry is not valid for high 
Tc superconductors Q , which in fact should contain all 
components of L = 0,2,4, .... Another argument comes 
from the newly discovered iron-based superconductors, 
which have a gap symmetry ~ cos k x cos k y [3, Q . This 
symmetry cannot be explained by the above s, d, g-wave 
interpretation. It is called an extended s-wave. Some re- 



*Electronic address: wliming@scnu.edu.cn 



searchers such as Yao et al claimed a d xy gap symmetry 
on one hand and at the same work they have to admit a 
sign change for the gap on T and jVfJjj in the Brillouin 
zone. It is seen that the s- or d-wave description becomes 
very controversial for the gap symmetries in iron-based 
superconductors. 

In this paper I argue that the gap symmetries of square 
lattice superconductors in fact originate from the 4-fold 
rotational symmetries of the coupling potential on neigh- 
bor sites. It has nothing to do with the angular momen- 
tum. I present a theoretical analysis and a computational 
demonstration for this argument. It is shown that the at- 
traction between the nearest neighbors (n.n.) or the next 
nearest neighbors (n.n.n.) provides a coupling for Cooper 
pairs, so that the gap functions of Cooper pairs take the 
symmetries of the coupling potential between neighbors. 
It is interesting to notice that most gap symmetries of 
square lattice superconductors can be reproduced in this 
mechanism. 



II. GAP SYMMETRIES 

The BCS coupling Hamiltonian is written as 

Hbcs = -V^24l c -kl c -k'lCk>t (1) 
hh 1 

Change it into the site configuration then one finds 

Hbcs = -V ^ 4Ai c oi c ^ ( 2 ) 



It comes from an on-site coupling, V(ri,r2) = V8{t\ — 
r 2)[9(- This coupling indicates that in fact a Cooper pair- 
ing can only occur on the same site. The electrons in 
cuprate superconductors, however, have strong on-site 
repulsion due to the strong correlation as seen in the 
Hubbard model. One usually rules out the double occu- 
pancy through the method of Gutzwiller preiection[Icil|. 
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Because of this strong correlation effect the conventional 
BCS on-site pairing is certainly not energetically prefer- 
ential in the cuprate superconductors. Hence, the inter- 
action Hamiltonian has to be changed to the coupling 
between neighbor sites as follows 

Hint = — ^ VmCi+mt C il C jl C j+m1 ( 3 ) 
ijm 

where V m is the coupling potential between two electrons 
(or holes) on two neighbor sites, here m denote the neigh- 
bors of a fixed site. In general, only the potential values 
on the near neighbors, such as Vq, V±, V2, are important. 
In the mean field approximation the above Hamiltonian 
becomes 

H mt = - ^(A£c- fei c feT + Afe4 T cL fci ) + E (4) 

k 

where A k is the gap function giving by 

A k = ~j2 v ™ A ™ e ~ ik ' Rm ( 5 ) 

771 
771 

With A m =< J2i c il c i+rnl >■ 

An important feature is that the range of the coupling 
potential gives the symmetries of the gap function, as 
listed in the Table I. 



Table I. Gap symmetries. 



m 


condition 


A k 


symmetry 







~ 1 


s 


n.n. 


(1) 


~ (cosk x + cos k y ) 


s~ 


n.n. 


(2) 


~ (cosk x — cos k y ) 


d x 2 — y 2 


n.n.n. 


(1) 


~ cos k x cos k y 


^x 2 +y 2 


n.n.n. 


(2) 


~ sin k x sin k y 


d X y 



In Table I n.n. denotes the nearest neighbor sites, 
x = (±1,0), y = (0, ±1), condition (1) stands for the 
isotropic case A x = A y , and (2) for the anisotropic case 
A x = —A y . In the case of the next nearest neighbor sites 
(denoted by n.n.n.), xy — (±l,±l),xy = (±1, =pl), con- 
dition (1) stands for the isotropic case A xy — A xy -, and 
(2) for the anisotropic case A xy = —A X y. In the above 
table the gap symmetry s x 2 +y 2 has been found recently 
in the iron-based superconductors @, Q]. Surprisingly, it 
is seen from the table that this neighbor coupling gives 
most symmetries of gap functions that have been found 
in various superconductors, including conventional and 
high-Tc superconductors with square lattices. 

It should be emphasized that these gap symmetries 
have nothing to do with angular momenta of the Cooper 
pairs. They originate from the point group symmetries of 
the crystals of superconductors. In fact, strictly speak- 
ing, the angular momentum L of Cooper pairs is not a 
conserved quantity since the potential does not obey the 
rotational symmetries of group i?3 0j which is only true 
for the local atomic orbitals. For the neighbor coupling 



the point group symmetries are much more significant 
than the atomic orbital symmetries. 

In order to confirm the isotropy and anisotropy of A& 
Green's function is defined as follows, 

« CitWMO) » = -i6(t) < [cn(t),c n (p)]+ > (?) 

According to the standard method of Green's function 
theory one obtains 

« c k] \c_ ki » u = & = J e 2 k + |A fc p (8) 

where e k are the energies of free electrons, e k — 
—2t(cosk x + cos k y ) + It' cos k x cos k y . Then the spec- 
trum theorem of Green's function gives 

< C-fcJ.Cfe T > = Im « Cfc T |c_ fei »u+iq 

7T 

tanh(q/2fc 5 r) 
= ^7 A k (9) 

2? k 

Combining the above equation and (O one obtains the 
gap equation 

A fc = -Enfe- g ) tanh( t /2fcBT) A 9 (10) 

where V{q) are the fourier components given by V(q) = 

( l / N )J2 m V ^ iq ' Rm - Under the 4 " fold rotations, R, of 
the point group symmetries of the potential U(m) this 
equation transforms in the following way (see Appendix), 

A Rk = -^V(k- q ) t&nH ^ 2kBT) A Rq (11) 

Compare the above equation with the gap equation pop 
one finds that solutions with the following two symme- 
tries exist: 

A Rk = ±A k (12) 

Therefore, both conditions (1) and (2) in Table I can 
be realized, that is, both the isotropic case and the 
anisotropic one are allowed. It is then confirmed that 
various gap symmetries listed in Table I are allowed by 
the neighbor coupling. A real system selects one of the 
symmetries with the lowest energy, or the largest gap 
magnitude. 

The symmetries revealed in (fT2"|) are called respec- 
tively conventional and unconventional superconductiv- 
ity in literatures 11]. In this paper it is proved that both 
these two symmetries are in fact a conclusion of the 4- 
fold rotational symmetries of the coupling potential on 
the square lattice superconductors. 

The gap symmetries can be checked by numerical 
calculations. The gap equation (fT0|) is solved self- 
consistently through iterations from a guessed initial gap 
function. It is found that this equation has stable solu- 
tions for a wide range of parameters. Two examples are 



3 



\ 




/ / « 










•Of— 










H V 

01 


° 


? 

( 








a-' t, y 




, o 






V \ 


-0.5 




k/]i 


0.5 



J i 


9 ° ? 

J g 


° V 


V, 
>J 




°'°> 0.01- *' 


v: 




\ \ 


O.Oj, 

/ r~~\\ 




-0.5 0.5 

k/x 



FIG. 1: Gap functions with the d x 2_ y 2 symmetry due to the 
n.n. coupling ti = 1.0, t 2 = 0.0, Uo = 0.0, f/i = -l.,U 2 = 
0.0, /i = 0.0(left panel) and with the s x 2 y 2 symmetry due to 
the n.n.n. coupling t\ = 1.0, ti = 0.0, Uo = 0.0, U\ = 0., Us = 
— 1.0, fj, = 0.0(right panel). 



shown in Fig.l for the d x 2_ y 2 and the s x 2 y 2 gap sym- 
metries. With the same parameters one also finds s~ 
and d xy gap symmetries from a different guessed initial 
gap function, but their magnitudes of gap functions are 
a few orders smaller. That is to say, the two symmetries 
shown in Fig.l arc most energetic preferential for the 
given parameters. Further computations show that dif- 
ferent guessed initial gap functions lead to different gap 
symmetries, but the gap function with the largest mag- 
nitude is almost determined by the neighbor coupling 
with lowest potential value. Thus the earlier theoretical 
analysis of this paper is further demonstrated by these 
calculations. It is concluded that the gap symmetries of 
cuprate superconductivity are in fact due to the neighbor 
coupling with 4-fold rotational symmetries of crystals of 
superconductors. 

The gap symmetries provide information for the cou- 
pling range of Cooper pairs. Comparing with the on- 
site coupling of the conventional s-wave BCS supercon- 
ductivity the Cooper pairs of high-Tc superconductors 
can distribute on two neighbors, including the n.n. and 
the n.n.n. neighbors. The size of the Cooper pairs is 
determined by the minimum of the coupling potential. 
The neighbor coupling probably provides a unified phe- 
nomenological model for the BCS and the high-Tc super- 
conductivity. 



III. CONCLUSION 



In conclusion, it is found that the gap symmetries are 
simply determined by the 4-fold rotational symmetries 
of the coupling potential on neighbor sites. A local on- 
site coupling potential results in the on-site pairing with 
the conventional s-wave symmetry, but a n.n. or n.n.n. 
coupling potential results in the pairing on neighbor sites 
with the s~ , d x 2_ y 2, s x 2 y 2 or d xy symmetries. It is proved 
that both isotropic and anisotropic gap functions are al- 
lowed by the 4-fold rotational symmetries of the coupling 
potential. Finally a numerical computation is performed 
to demonstrate the gap symmetries. The solutions of the 
gap equation are very stable for a wide range of param- 
eters. 
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Appendix 

Eq. (jlip can be proved in the following way. From 
(fTUll one has 



tanh(^/2fc B T) 



A Rk = 



qm 



,e" v " H " — -A„ 

2( 

N 2_, VRm€ 2 £ 9 



qm 



\ ^ Tr i(q — k)-m 

tanh(^/2fc B T) 



Rq 



Rq 



where the symmetry m £ g has been used since the 
value of ^ q depends weakly on A q and en q = e q . 



[1] Gabriel Kotliar and Jialin Liu, Phys. Rev. B 38, 5142 
(1988) 

[2] C.C. Tsui et al. Nature 387, 481 (97) 

[3] Yunping Wang, http://www.sciencenet.en/m/user_content 

aspx?id=211798. 
[4] Yan-Yang Zhang, et al, Phys. Rev. B 80, 094528(2009). 
[5] M.M. Parish, et al, Phys. Rev. Lett. 78, 144514 (2008). 
[6] LI. Mazin, D.J. Singh, M.D. Johannes, and M.H. Du, 

Phys. Rev. Lett. 101, 057003(2008). 

[7] zi-Jian Yao, Jian-xin Li, Z.D. Wang, larXiv:0804.4l66V 3 



[cond- mat. supr-con] 2008. 
[8] S. L. Yu, J. Kang, and J. X. Li, Phys. Rev. B 79, 064517 
(2009). 

[9] J. F. Annett, in Superconductivity, Superfluids and Con- 
densates(ed. by Oxford Univ. Press 2004),pl34. 
[10] Attaccalite Claudio, Fabrizio Michele, Phys. Rev. B68, 
1551171(2003) 

[11] J. F. Annett, in Superconductivity, Superfluids and Con- 
densates(ed. by Oxford Univ. Press 2004),pl58. 



